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A density dependent, effective nucleon-nucleon force of the Skyrme type is derived from the quark- 
' meson coupling model - a self-consistent, relativistic quark level description of nuclear matter. This 

new formulation requires no assumption that the mean scalar field is small and hence constitutes 
, a significant advance over earlier work. The similarity of the effective interaction to the widely 

£\j ■ used SkM* force encourages us to apply it to a wide range of nuclear problems, beginning with the 

binding energies and charge distributions of doubly magic nuclei. Finding acceptable results in this 
conventional arena, we apply the same effective interaction, within the Hartree-Fock-Bogoliubov 
approach, to the properties of nuclei far from stability. The resulting two neutron drip lines and 
shell quenching are quite satisfactory. Finally, we apply the relativistic formulation to the properties 
of dense nuclear matter in anticipation of future application to the properties of neutron stars. 
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I. INTRODUCTION 

o 

m 

Given that QCD is almost universally believed to be the fundamental theory of the strong interaction, one of the 
\Q , great challenges facing modern nuclear physics is to derive the properties of atomic nuclei from it. Amongst the 
hopes for a more fundamental description of nuclear matter is to be able to extrapolate to extreme regions of density, 
asymmetry or even strangeness with rather more confidence than has hitherto been possible. One might also expect 
to understand the nuclear EMC effect 0, 0, and a number of other challenges to traditional nuclear theory 

While at the present time there has only been exploratory work aimed at relating QCD itself to nuclear structure _ 
0] , there has been considerably more work within the framework of quark models and effective field theories utilizing 
the symmetries of QCD. Here we shall concentrate on one particular model, the Quark-Meson Coupling (QMC) 
model 0, 0, El , which is built around the self-consistent response of the relativistic, confined valence quarks to 
the Lorentz scalar mean field in nuclear matter. In previous work we have shown that the main features of the 
effective interactions used in nuclear physics can be quantitatively derived from the QMC model. In order to do 
this we rewrote the QMC model as an equivalent, effective Hamiltonian with zero range many body interactions. 
The many body character of the interactions was a direct consequence of the response of the quark structure to the 
nuclear environment. We found that this Hamiltonian was in quite good agreement with the popular Skyrme-III 
interaction [12|, which led us to conclude that quark level dynamics does play a key role in ordinary nuclear physics, 
even if it can be hidden within an effective many-body theory. 

However, one unsatisfactory feature of the derivation in Ref. [Tl| was that the approximations which were required 
meant that one could not justify the use of the resulting effective Hamiltonian much beyond the density of ordinary 
nuclear matter - (po — 0.16fm~ 3 ). One of those approximations is actually tied to the physics of the model: the 
response of the quark structure to the nuclear medium leads to a non-linear meson field equation which is intractable 
in its full generality. In Ref. |Tlj we relied on an iterative solution of the field equations, which essentially amounted 
to an expansion in powers of the density. We showed that this is an acceptable approximation for ordinary nuclear 
densities, with the medium effects then appearing as many body forces. Another approximation used in that work 
was a non-relativistic expansion which, in the framework of the iterative solution, appeared unavoidable. The latter 
approximation was also motivate d by our wish to compare the effective Hamiltonian with the Skyrme interaction 
which is a non-relativistic object. |26| These approximations, which were justified in the context of the earlier work, 
become untenable at higher densities. For instance, in this model one finds that the velocity of sound exceeds the 
velocity of light at about 3po- This is somewhat frustrating because it prevents us applying the model to the study of 
neutron stars. 

In this work we propose a new formulation which does not rely on the approximations which we have just outlined. 
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Instead of solving the meson field equation iteratively, we expand the field around its nuclear expectation value and 
treat the deviation as a perturbation. This leads to an effective Hamiltonian with a density dependent two-body 
force. As the equation for the expectation value can be solved exactly (at least in a numerical sense), the bulk of the 
medium effects may be taken into account without any assumption about the value of the density. This is a major 
advance with respect to the formulation in Ref. Moreover, the case of uniform nuclear matter, which is the 

relevant approximation for neutron stars, can now be treated in a way consistent with relativity. 

This paper is organized as follows. In Section[H]we recall the basic elements of the QMC model. This is followed in 
Section fllll bv the formal derivation of the equivalent, density dependent NN force. A non-relativistic reduction of the 
full effective interaction is carried out in Section lTvl in order to make a detailed comparison with the density dependent 
Skyrme force in ScctiorlVI In Section lVTl we carry out a full Hartree-Fock calculation with the non-relativistic, effective 
force that we derived and compare the results for selected finite nuclei against experimental data. Finally, in Section 
IVIII we solve the full model for uniform matter containing just nucleons - leaving for future work the generalization 
of the method to include hyperons. Section [Villi presents some concluding remarks and suggestions for further work. 



In the Quark Meson Coupling model, the nuclear system is represented as a collection of confined clusters of three 
valence quarks. Because it is possible to carry through much of the calculation analytically, the MIT bag model was 
chosen to confine the quarks. It is assumed that the effect of having these bags overlap, that is clusters of more than 
three quarks, can be neglected. This approximation certainly breaks down at some density but we assume that it is 
still acceptable in the density range that we consider. In this respect, we observe that the bag model is an effective 
realisation of confinement, which must not be taken too literally. Indeed, QCD lattice simulations 0, 0] strongly 
suggest that the true picture is closer to a Y-shaped color string attached to the quarks. Outside this relatively 
thin string one has the ordinary, non-perturbative QCD vacuum where the quarks from the other nucleons can pass 
without disturbing the structure. Thus, while the bag model imposes a strict boundary condition which prevents 
the quarks from travelling through its boundary, this must be seen as the average representation of a more complex 
situation. One should not attribute a deep physical meaning to the boundary of the cavity nor to its size. In particular, 
estimating the density at which the non-overlap approximation breaks down as p — l/V, with V the bag volume, is 
certainly too pessimistic. 

The salient feature of the QMC model is that the interactions are generated by the exchange of mesons coupled 
locally to the quarks. In a literal interpretation of the bag model, where only quarks and gluons can live inside the 
cavity, this coupling would be unnatural. On the other hand, in terms of the more realistic picture presented by lattice 
QCD, the quarks are attached to a string but otherwise move in the non-perturbative QCD vacuum. There nothing 
prevents them from feeling the vacuum fluctuations which we describe by meson fields. As in our earlier work, we 
limit our considerations to the er, uo and p mesons. We recall that the a meson here is not the chiral partner of the 
pion. It is a chiral invariant field which represents correlated two-pion exchange. The quarks of the bag are in the 
Weinberg representation |l5l llq , where the mass term qq is a chiral scalar which allows a chiral invariant coupling of 
the form aqq. 

Of course, this is not the end of the story because the transformation to the Weinberg representation introduces 
an infinite set of derivative pion-quark couplings. These produce a complicated set of multi-pion exchanges, together 
with contact interactions. We assume, in the spirit of both Quantum Hadrodynamics 17] and the one boson exchange 
model, that the bulk of the corresponding physics can be modeled by the cr, u> and p exchanges between quarks. It is 
clear that the long ranged single pion exchange does not fit in this picture and must be treated separately. However, 
because of its pseudoscalar nature, this exchange does not contribute to the mean nuclear field. In the Hartree Fock 
approximation , which will be our framework in the following, it contributes through its exchange (Fock) term. This 
effect has been estimated in the framework of infinite nuclear matter and found to be small. More exactly, most of 
the effect can be absorbed by a readjustment of the couplings to the other mesons. This is our first motivation for 
neglecting this part of the interaction but we must keep in mind that the argument is only true in the Hartree Fock 
approximation. When other nucleon correlations are taken into account, the game changes but this is beyond the 
scope of this work. 

Our starting point is the classical energy of a system of non-overlapping bags of quarks, which are coupled to the 
nuclear meson fields a, lo and p. Following Ref. |ll| we can write 



where the first term is the relativistic energy of nucleon i, with momentum Pi and effective mass M e tf, g^u is the 
repulsion felt by each nucleon as a result of its vector interaction with the time component of the lo mean field, V so 



II. THE QMC MODEL 
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is the spin-orbit potential and the static meson energy is 

dr 



with m a ,m u ,m p the meson masses. As usual, we consider only the time component of the vector fields. The effect 
of the p, which can be treated by analogy with the u> field, will be introduced at the end. The expression for the spin 
orbit interaction, V so , first derived in Ref. 0, is not necessary at this stage. 

The effective mass, M e ff(a), is the rest frame energy of a quark bag in the a field of the medium, evaluated at the 
center of the bag. In order to calculate M e ff(<r) one needs to solve the bag equations for the relevant value of the a 
field. However, we have checked that a quadratic expansion: 

M eff (a) =M~g a a+^(g a a) 2 (3) 

is accurate up to values of g a a as large as 600MeV, which corresponds to densities far beyond the physical limitations 
of the model. The parameter d is the scalar polarizability of the nucleon. It depends explicitly on the response of the 
quark structure to the external scalar field. In the bag model it is well represented as a function of the bag radius, 
Rb, as 

d= 0.0044 + 0.21LR B -0.0357i?|, (4) 

where both d and Rb are in fm. This does not give exactly the values that were used in Ref. because we have 
now included the contribution of the spin dependent, "hyperfine" color interaction to the bag energy. The coupling 
constant, g a , which by definition refers to the nucleon, is related to the cr-quark coupling, g%, by 

9<j = 3s£ / drqq{r) , (5) 

J Bag 

where q is the valence quark wavefunction for the free bag. For the vector couplings the relationship is 

9u = 3gl, g P = g q p . (6) 

Equation @ holds for any flavor, /, of hadrons, provided the couplings and scalar polarizability are replaced by their 
corresponding values g<r(f),d(f). In the bag model they are related to g a and d in a well defined way. Thus the 
treatment of flavors other than the nucleon is straightforward and does not introduce new parameters. This matter 
will be developed in forthcoming work. Here we limit our consideration to the case where the medium contains just 
clusters with the quantum numbers of protons and neutrons. 

III. HAMILTONIAN 

By hypothesis, the meson fields are time independent. Therefore the classical Hamiltonian for the nuclear system 
is simply 

H(Ri,Pi) = E(Ri, Pi,a ^ a soU u ^ u sol ) , (7) 

where Rj and Pi are the position and momentum of nucleon i and a so i , lo so i are the the solutions of the equations of 
motion corresponding to Eq. l|T|l 

SE _ SE _ 
So- Slu 

that is 

-VV + m> = -Y^Sif-Rt^Pl + Meffiay, (8) 

i 

-V 2 Lo + mliu = g^Sp-Mi). (9) 

Note that in Eqs. (|8I9|I we have neglected the contribution corresponding to the variation of the spin orbit interaction 
V so - As pointed out in Ref. this results in an error of order Vf in the Hamiltonian and it is consistent to neglect 
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it because the spin-orbit interaction in Eq. has been derived as a first order perturbation Since the equation 
for the D field is linear, its solution is elementary and poses no new problem with respect to our previous work. Its 
contribution to the energy will be written later and we now concentrate on the field equation for the a. 
We assume that it makes sense to write 

a =< a > +5a , 

where the C-number < a >, also written a, denotes the nuclear ground state expectation value, that is p7) 



< a(r) >= J dR 1 ...dR A <f>*(R 1 ...R A )a(r,R l ,P l )<f>(R 1 ...R A ) 
and to consider the deviation, 5a, as a small quantity. If we define 



K = J2S(r - R t )^/P 2 + M eff (a) 2 , 

i 

we see that the a field equation has the form 

/ -2 2\ / e n 9K dK , ^ . d 2 K . . 
(-V 2 + ml) (a + Sa) = -— = - — (a) Sa^(a) 



da da 

We also expand ^ (a) , ^Jf (a) about their expectation values 

dK. , dK, . 

— (a) = < —(a) > +5 

d 2 K. . d 2 K . . 



and suppose that 



<5<r, 5 



dK , 



d 2 K 
da 2 



d 2 K 
da 2 



(<?) 



are small quantities. We can then solve the meson field equation order by order, which gives: 

(-v 2 



'm 2 ) a 



dK,. 



-V 2 + ml) 6a 



-6 



dK 
~da~ 



6a < -—{ a )>, 
da z 



dK 
~~da~ 



dK 



d 2 K 



(a)+ < > -6a < -5-5- (o-) > 



0. 



dc 



(10) 



(11) 



(12) 



(13) 



(14) 



As we limit the expansion of the Hamiltonian to order (5a) 2 , it is sufficient to solve the field equation at order 6a, 
which corresponds to Eqs. i|13ll4|) . Using integration by parts, the Hamiltonian may be expanded as 



H = J d rK\ s +5a^(a)+ 1 -(5a) 2 ^(a) 



+ a (-V 2 + m 2 a ) a + 5a (-V 2 + m 2 ) a + ]-5a (-V 2 + m 2 ) 5a . 

— _ 



(15) 



To this order we can replace 



d 2 K._. d 2 K.. 



da 2 



dc 



because this multiplies (5a) 2 . Using Eqs. i|13ll4fl we then find: 



H = dr 



1 



dK 



1 



K(a) --*< —(a) > +-6a —(a)- < —(a) > 



da 



dK 



dK 



da 



da 



(16) 
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Note that the mean field approximation amounts to neglecting 8a in Ea. (|16|l . To complete the derivation we need a 
prescription for writing the quantum form of K and its derivatives. The important simplification is that these are 
one body operators because they are evaluated at the C-number point a = a. Thus we can write 

a/3 

where a^ a ,a a are the creation and destruction operators for the complete 1-body basis \a >. The matrix elements 
K a p(a) must be chosen so as to reproduce the classical limit, Ea. (|l(JII . In the momentum space representation, there 
is a natural choice 

K & = h 5 eiM3 ' ? (v /fc2 + Me ^^^ 2 + + M eff [a(r)}^ a\a % , , (17) 

k,k' 

where the symmetrization is introduced to guaranty hemiticity and V is the normalisation volume. We also choose 

^^^''l (^fc 2 + M e}f [a{r)Y + V^' 2 + M eff [a(r)}^ a^p , (18) 

k,k' 

with a similar expression for the second derivative. The ordering ambiguities associated with products of non- 
commuting operators are fixed by the normal ordering prescription, which amounts to removing that part of the 
energy which originates from the interaction of one nucleon with its own field. 

The Hamiltoniam defined in Eq. I|ltj|) is not a standard many-body problem, since we do not know a and 8a until 
the ground state nuclear wave function has been specified. Therefore the Hamiltonian must be determined (through 
a and 8a ) at each step of the self consistent procedure. This is a significant technical complication. On the other 
hand, for our purposes it is not necessary to solve this Hamiltonian in the general case. One of our main goals is 
to obtain the equation of state for very dense nuclear matter, with a view to applying it to neutron stars |18| . but 
in this case we only need the uniform matter approximation for which the model is easy to solve. Of course, we do 
need to consider finite size effects for ordinary nuclei, a necessary step in order to show that the model is realistic 
with respect to nuclear phenomenology. However, in this case, we can make a non-relativistic expansion and build a 
density functional, < H{f) >, which can then be used in variational calculations of nuclei. 



IV. NON RELATIVISTIC EXPANSION 



We first consider the case of finite nuclei, in order to verify that the model is indeed realistic. For this purpose we 
build the density functional, < H(r) >, using approximations which are standard in low energy nuclear physics. The 
approximations that we use below involve a non-relativistic expansion and the neglect of those velocity dependent 
and finite range forces which involve more than 2 bodies - as is the case for conventional effective nuclear forces. This 
amounts to expanding the terms which involve either the momentum or the gradient of the density in powers of the 
a nucleon coupling and stopping at order g 2 . To simplify the expressions we omit the spin and flavor indices as long 
as they are not truly necessary. We first define the number density, D(r), and kinetic density, £(r), by 



D{r) 



1 

V 



E 



k.k' 



e Kk-k').f a l 

k K 



i 



E i(k-S').f k 
k,k' 



(19) 



Then, following the approximation scheme defined above and using Eq. © we find the following expression for the 
operator K and its derivatives 



K\, = D(r)M e// [^] + ^(l + ^), 



dK 
~da~ 

d 2 K 



BM, 



8a 2 



D{f) 



= dglD(r) 



eff 



9c 



m 

2M 2 ' 



(20) 



Substituting these expression into Eqs. (|13I14|) and using the same approximations, we can solve for a and 8a: 



9*° 



G a < D > 
1 + dG a < D > 



G, 



2M 2 



G(j 



V 2 < D > 



(21) 
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5a = 



1 / dK 



(*)+<^(*)>)+^v^ 



1 „ 2 1 / dK, . dK . . . 

^(a)+<^(a)>), 



9cr No- 
where we have defined G CT = g^/m 2 and the (position dependent) effective a mass 

m 2 , = ml (1 + Gad <D>) . 



(22) 



(23) 



We can now substitute Eqs. I|2UI21I22|I into Eq. nitty to obtain the Hamiltonian of the model. In practice we want the 
corresponding density functional in the Hartree Fock approximation, so we need to evaluate (note that by definition 
< 5a >= 0) 



TT,-* 1 dK . . 1 . dK,. 

< H(r) > = < K(a) > --a < —(a) > +- < 5a— (a) > . 



(24) 



where the ground state wave function is a Slater determinant, with Fermi level F, built from the single particle wave 
functions $ J (r, a, m). We now restore the spin flavor dependence but restrict our considerations to nuclei made of 
protons and neutrons. So the flavor index is just the isospin projection m = ±1/2. As usual, we define |l2j | : 



Pm (r) = E E I $ 4 (r , cr, m) \ 2 , p(r) = E p m (r) , 

T ™ if) = E E I ^ TO °) 2 ' T ^ ) = E T ™ (*0 = 



«GF cr 

Jm(f) = i EE^'' 7 

iGF crcr' 

and using standard techniques we find, after some algebra: 
<H(f)> = pM+^ 



V$ l (r, a, m)] (r, cr', m) , J{r) = E -4 ( f ) ■ 



(25) 
(26) 
(27) 



^2 (p t +^E^" v2 ^) 

pV 2 p - E ( JPm^Pm - PmT m J J 
m / 



2M 2 

I G a G a 
1 G, 



[f 2{l + dG a p)^ pl ) 



2 1 + dG a p 

For completeness, we give the ui,p and spin orbit contributions (labelled IS for isoscalar and IV for isovector): 



<: H*J{r) :> 
<:Hfr^:> 
<: W w (f) :> 



1 



4M 2 
4M 2 



G a + ( 2—— — 1 I G u 



/<iV 



E 



1 + ^mm' ] Pm'V.Jn 



2^-1 



E 



mm' + -C mm i ) p m '^-Jn 



(28) 

(29) 
(30) 



G' 



E 



1 2 ^mm' ) PmPri 



2 ( P m ^ Pm' , &mm' iyPm^ Pm ^Pm^~rnj 



<: H p {r) :> = — ^ E j f™' _ \c mm , J p m/ o m ' 



(31) 



(32) 

//(,- v * J ">-~ P ) 

where: C mm / = S mm im 2 + {5 m , m '+\ + 5m',m+i)/2 and G^ = g 2 / m 5> ^p = 9p/ m % The isoscalar and isovector magnetic 
moments which appear in the spin orbit interaction have the values 

Pis = P P + Pn = 0.88, p IV = p P - p n = 4.7 . (33) 



„ | ill 772 A Cf nrn ' J PmV Pm' H~ o Ctnrn' PmTn 
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771a (MeV) 


G CT (fm 2 ) 


G^fm 2 ) 


G p (fm 2 ) 


K N {MeV) 


600 


12.652 


9.838 


9.67 


346 


650 


12.428 


9.308 


8.583 


346 


700 


12.254 


8.899 


7.724 


346 


750 


12.116 


8.575 


7.048 


346 



TABLE I: The couplings GcG^yGp for a bag radius Rb = 0.8fm and several values of the scalar meson mass, m a . The last 
column is the nuclear incompressibility. 



As we have already pointed out, these expressions are the same as in our previous work. We note that all terms which 
involve the square of the spin density, J, have been neglected. This is common practice and, since it amounts to 
treating the spin orbit interaction as a first order perturbation, it is consistent with our derivation of the expression 
for the effective mass, Eq. (Q. 

It is convenient to have a more explicit expression for the density functional. Using Eqs. 1)281- 1^1 one can write, 
using the same notation as in Ref. 



cxpli 

m 



< H(f) >-- 



PM+ 2M 



Ti-o + Hz + H e ff + Hfi 



(34) 



where 



7~to + W3 



n 



eff 



Hfi 



~3G P j 

32 8(1- 

5 Gp 



(Pn ~ Pp) 



Gfj G & 3 G^ 

hd P G a ) 3 ~ 2 (l + d P G a ) ~~ 



G„ 



G u 



32 8(1 + d P G a f 8 



G a 



G ( 



\^8m p 2 2m a 2 2m J 2 
-p <-> n, 



Go \ ( G p 

4M N *J Pn+ Um p z 



G, 



m?) PP . 



r„ 



3 Gp 3 Go 3 G^ G { 



y32m p 2 8m a 2 



mj 1 



3 Gp G, 



v 16m p 2 2 



(7 . G^ 
ma- 2 



In 
G a 



V- J n 



2 m J AM N 2 ) Pp . 
3 G a 3 G^ (—1 



-G a 



1M N 



1 .\/\ :! i .\/.\ 



n 2 8 Mn 2 

Gcj (l-2p s ) \ 
a K/r. 2 p P 



In . 
Ips) _ 3_G, 



V 2 (p„) +p <-> n, 



2 (Pn 

-1 + 2 



32 M N 2 



Pn 



+ p <-> 77. 



(35) 
(36) 

(37) 
(38) 

(39) 

(40) 

(41) 



We have determined the couplings G a , G^, , G p by fixing the saturation density and binding energy of normal nuclear 
matter to be po = 0.16fm 3 and Eb = —15.85 MeV, as well as the asymmetry energy of nuclear matter as 0,4 = 30 
MeV. Apart from a small readjustment of the couplings, we found no significant sensitivity to the bag radius. We 
therefore display our results for only one value, Rb = 0.8fm, which is quite realistic ^(|. The u> and p masses are 
set at their experimental values. The last parameter, which is not well fixed by experiment, is the a mass. We shall 
use m a — 600, 650, 700, 750MeV. The corresponding results are given in Table U We see that the incompressibility, 
K N , is a little high with respect to the currently prefered range, but we point out that this calculation has not yet 
taken into account the single pion exchange interaction. We know that the pion Fock term alone reduces Kn by 
as much as 10% and it is likely that this is amplified by other correlations. 



V. COMPARISON WITH THE SKYRME FORCE 



As a means of orientation, we compare the density functional of our model with that corresponding to a typical 
density-dependent effective interaction. We choose to compare with the popular Skyrme SkM* parametrisation rather 
than more recent ones, such Sly4 0, because in the latter case the formal identification is complicated by the larger 
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m CT (MeV) 


io(fm 2 ) 


ti(fm 4 ) 


t 2 (fm 4 ) 


Him 5 / 2 ) 


xo 


W)(fm 4 ) 


Deviation 


600 


-12.72 


2.64 


-1.12 


74.25 


0.17 


0.6 


33% 


650 


-12.48 


2.21 


-0.77 


71.73 


0.13 


0.56 


18% 


700 


-12.31 


1.88 


-0.49 


69.8 


0.1 


0.53 


18% 


750 


-12.18 


1.62 


-0.28 


68.28 


0.08 


0.51 


38% 


SkM* 


-13.4 


2.08 


-0.68 


79 


0.09 


0.66 


0% 



TABLE II: Comparison of the SkM* parameters with the QMC predictions for several values of m a 



number of parameters. The SkM* interaction depends on 6 parameters: to, t\ 1 t 2 , £3, xq and Wo and its energy 
density |19| | may be written, using the same notation as in the previous section: 



Ho + H3 
n eff 

Hfi n 
H so 



pit 3 (2 p 2 - Pn 2 - Pp 2 ) t (p 2 (2 + x ) - (1 + 2 s ) (pn 2 +P P 2 )) 

24 + 4 

pr (2ti + 2t 2 ) (-h + t 2 ) {p n r n + p p t p ) 

8 8 
-(p (6tx-2t 2 ) V 2 (p)) (3h+t 2 ) (p„V 2 ( l o»)+P P V 2 (Pp)) 



32 

(Wo GoV ■• J + V • J„p„ + V • J p p p )) 



32 



(42) 
(43) 
(44) 
(45) 



To simplify, we compare H e ff,Hfi n and H so with the QMC expressions Eqs. (|37I39I41(I . in the case N = Z. This 
detemines tt,ts and Wo- To find to,t3 and Xq, we consider the term H0 + H3. As the functional form is not the same, 
we fit (Ho + Hz)q MC with the form which appears in (Ho + Hz) skM in the range p E [0 — > 0.2 fin -3 ]. We first do 
this in the case N = Z, which determines to and t$ and using these values we do the fit for Z/A = 92/208, which 
determines xo- The results are collected in Table ITT1 for several values of the a mass. We find satisfactory agreement 
with the SkM* parameters in the window m a = 650 + 700 MeV. This comparison suggests that our model may provide 
an acceptable representation of low energy nuclear physics. 

However, this comparison is really rather qualitative. First, we have minimized the isospin effects by setting N = Z 
when determining the parameters t\,t 2 ,Wo- Second, the values we find for to and t% depend somewhat on the range 
of density we use for the fit. Thus a more direct comparison of our model with actual nuclear data is desirable. This 
is done in the next section. 



VI. HARTREE FOCK CALCULATIONS FOR FINITE NUCLEI 



As pointed out in the previous section, the QMC and the Skyrme energy functionals have a similar structure. 
However, the density and the isospin dependence of some terms, particularly H3 and H so , are rather different in the 
two approaches. Thus, for the Skyrme functionals the term H3 has a density dependence of the form p a . Originally 
this term was interpreted as having been generated by a three-body contact interaction, which is equivalent to a 
two-body density-dependent force in even-even nuclei [lj] ■ In the QMC functional, the density dependence of H 3 is 
much more complicated. However, a more complicated density dependence could also appears in Skyrme type energy 
functionals if they are derived from a more fundamental theory based on G-matrix, such as density-matrix expansion 
(DME) method of Negele and Vautherin |20|. A comparison between the energy functionals corresponding to DME 
and to standard Skyrme forces shows that in fact the strength of the density-dependent term is rather similar in 
both cases, in spite of their different density dependence. It is interesting to observe that a density dependence of a 
fractional type, as appears in QMC, was previously considered in phenomenological functionals by Fayans et al. |2l| . 
The advantage of fractional expressions in particle density is that the corresponding energy functional preserves causal 
behaviour up to very high densities. This is actually the case for the QMC functional as well. 

Another difference between the QMC and standard Skyrme energy functionals arises in the isospin dependence of 
the spin-orbit term H so . In both cases the form factor of the one-body spin-orbit interaction for a nucleon with the 
isospin m (n or p) can be written as: 

W m (r) = aVp m + bVpfh , (46) 
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where p m is the particle density, while fh denotes the opposite isospin to m. For the standard Skyrme forces the ratio 
s = a/b is equal to 2. This isospin dependence of the one-body spin-orbit potential is induced by the exchange term 
(since the Skyrme two-body spin-orbit force is isospin independent). For the QMC functional used in the calculations 
below, the ratio s is about 2.78. This strong isospin dependence in both the QMC and Skyrme functionals is in 
contrast with the weak isospin dependence employed in the relativistic mean field models, in which the contribution 
of Fock (exchange) terms is neglected • 

Starting from the QMC energy functional one can easily derive the corresponding Hartree-Fock (HF) equations. 
They have a form similar to the Skyrme-HF equations, apart from the rearrangement term and the one-body spin- 
orbit interaction, which (as discussed above) have a different densit y an d isospin dependence. The HF equations were 
solved in coordinate space, following the method described in Ref. [IJ and the Coulomb interaction was treated in a 
standard way - i.e., the contribution of its exchange part was calculated in the Slater approximation. The calculations 
were performed for the doubly magic nuclei 16 0, 40 Ca, 48 Ca and 208 Pb. For defmiteness, the a meson mass has been 
set to m CT =700MeV, as suggested by the comparison with the SkM* interaction. At this point we recall that the 
QMC model is essentially classical because both the position and velocity of the bag are assumed known in the energy 
expression . The quantization then leaves some arbitrariness in the ordering of the momentum dependent pieces 
of the interaction. As pointed out in previous work [Tl| . in the non-relativistic approximation the difference between 
the orderings is equivalent to a change of about 10% in m a . In this work the ordering is fixed by the relativistic 
expression chosen for the operator K, Eq. I|17l) . The non-relativistic reduction then leads to an ordering which is not 
the same as in Ref. [Tl|. This is why the a meson mass that we use here is somewhat higher. Note that this ordering 
ambiguity is only of concern in the case of finite nuclei. In uniform matter, which is the relevant approximation for 
neutron stars, the problem does not exist. 

The results for the binding energies and charge radii are shown in Table 3. The charge densities are calculated with 
the proton form factor usually employed in the Skyrme-HF calculations [l2( • From Table 3 one can see that QMC-HF 
gives results which are in reasonable agreement with the experimental values. The agreement is not as good as that 
given by the recent Skyrme or RMF models, but one should keep in mind that in these models the experimental values 
for the binding energies and radii are included in the fitting procedure, which is not the case for the QMC functional. 

One also finds a similarly reasonable description for the spin-orbit splittings, shown in Table 4. Since the isospin 
dependence in QMC-HF is stronger than in Skyrme-HF, one would expect different values for nuclei with large 
isospin asymetry. However, as can be seen from Table 4, the differences are rather small. This is primarily because 
the spin-orbit splitting depends on the product of the spin-orbit form factor and the corresponding single-particle 
wave functions. Thus, if the wave functions are not strongly localised in the surface region, where W T (r) is effective, 
the influence of the isospin dependence of W T (r) upon the splitting need not be so significant. 

In Figs, n 121 we show the proton and neutron densities calculated with the QMC model and with the Sly4 Skyrme 
force . In the proton case we also show the experimental values [23j . Once again the two models give rather similar 
results, with the largest differences noticed for the neutron skin of 208 Pb. More precisely, for this nucleus the neutron 
skin (i.e., the difference between the neutron and proton rms radii) is equal to 0.12fm for QMC, compared to O.I6fm 
for the Skyrme force. With respect to experiment the QMC model tends to overestimate the density in the central 
region, but the overall agreement is quite good, given that the model has no parameter adjusted to fit the properties 
of finite nuclei. 

Next we analyse the predictions of the QMC functional in nuclei far from stability. In order to do such a study, one 
must take into account the contribution of pairing correlations. They are treated here in the Hartree-Fock-Bogoliubov 
(HFB) approach. For the pairing interaction we have taken a density-dependent contact interaction of the form [24| : 



V = V 



V 



p_ 

Pi) 



6(n ~ r 2 ) ■ (47) 



With this pairing force the QMC- HFB equations are local. They have been solved in coordinate space by truncating 
the quasiparticle spectrum at an energy equal to 60 MeV. The parameters of the pairing force have been fixed to the 
values: Vq — —333 MeV fm -3 , rj = 0.5, a = 1, and p = O.I6fm~ 3 , which were chosen to yield a reasonable average 
gap for the tin isotopes ( a benchmark for pairing models). 

With this force we have investigated the position of the drip lines, which give the limits for bound nuclear systems. 
Here we present the results for the two-neutron drip line predicted by the QMC functional for Ni and Zr isotopes. The 
two-neutron drip line is defined by the value of N (number of neutrons) for which the two-neutron separation energy 
becomes negative. In the QMC-HFB calculations the neutron drip line appears around N=60 for Ni and around N=82 
for Zr. These positions of the drip line are similar to the predictions provided by the Skyrme force SLy4, commonly 
used in non-relativistic calculations [l9j|. 

Another important property which characterizes nuclei far from stability is the shell quenching, which has important 
consequences for the astrophysical rapid capture processes. The phenomenological energy functionals show that when 
one approaches the drip line the shell gaps associated with the magic numbers of neutrons or protons are quenched in 



10 





E B (MeV, exp) 


E B (MeV, QMC) 


r c (fm, exp) 


r c (fm, QMC) 


10 Q 


7.976 


7.618 


2.73 


2.702 


40 Ca 


8.551 


8.213 


3.485 


3.415 


48 Ca 


8.666 


8.343 


3.484 


3.468 


208 Pb 


7.867 


7.515 


5.5 


5.42 



TABLE III: Binding energy and radii calculated in QMC-HF, as described in the text. 





Neutrons (exp) 


Neutrons (QMC) 


Protons (exp) 


Protons (QMC) 


16 0, l Pl/2 - lp 3 /2 


6.10 


6.01 


6.3 


5.9 


40 Ca, lds/a - ld 5/2 


6.15 


6.41 


6.00 


6.24 


i8 Ca, ld 3/2 - lds/a 


6.05 (Sly4) 


5.64 


6.06 (Sly4) 


5.59 


208 Pfe, 2d 3 /2 " 2d 5 /2 


2.15 (Sly4) 


2.04 


1.87 (Sly4) 


1.74 



TABLE IV: Comparison between the QMC and "experimental" spin orbit splittings. Because the experimental splittings are 
no so well known in the case of 48 Co and 208 Pb, we give the values corresponding to the Skyrme Sly4 prediction. 



comparison with stable nuclei. The shell quenching can be tested by calculating the two-neutron separation energies, 
S% n , across a magic number. Here we present the QMC-HFB results for the neutron magic number N=28. The 
variation of S2n across N=28 is calculated for two extreme values of proton numbers, namely Z=32 (proton drip-line 
region) and Z=14 (neutron drip-line region). One thus finds that S2 n changes by about 8 MeV for Z=32 and by about 
2-3 MeV for Z=14. This strong shell quenching is very close to that obtained in the Skyrme-HFB calculations (see 
Fig. 15 of Ref. [13 ). 

In conclusion, the QMC functional gives a reasonable description of the known double magic nuclei and predicts, 
for nuclei far from stability, similar properties to those found in mean field Skyrme models. 



VII. HIGH DENSITY UNIFORM MATTER 



In the previous section we have shown that the non-relativistic approximation to the Hamiltonian developed in 
section Mil gives an acceptable phenomenology for ordinary nuclei. This give us some confidence to examine the 
consequences of the model in the high density region, relevant to the neutron star problem. The classical application 
of a density dependent force of the Skyrme type to this problem was that of Negele and Vautherin (25|. Of course, all 
non-relativistic effective theories must break down as the density increases and we therefore concentrate on a truly 
relativistic treatment in this case. For neutron stars one can restrict consideration to uniform nuclear matter, which 
implies that < K(a) > and its derivatives are numbers independent of r. The constant expectation value of the sigma 
field a is then determined by the self consistent equation (see Eq. JTSJl): 



1 dK,, 

— < > 
mi da 



(48) 



which is solved numerically. The fluctuation, da, is given by Eq. I|14|) . where the effective a mass, rh 
and the solution, expressed in terms of the a- meson Green function, Q a , is 



is now constant 



<5(T (r) 



dK 
~da~ 



dK 

< > 

OCT 



= / dr 



dq 



a iq.(r-r') 



(2tt) 3 q 2 + ml 



dK, dK, . 



(49) 



Note that, even though a is constant, the operator K and its derivative have an explicit dependence on position (see 
Eqs. (|17I18J| ) which disappears only in the expectation value. 
For uniform matter the energy density is 



e = = — I dr 

V V 



1 dK , . 1 . ,.dK. , 
< K(a) > --a < —(a) > +- < 5a(r) ^(^,r) > 



(50) 
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G CT (fm 2 ) 


G^fm 2 ) 


G p (fm 2 ) 


A"iv(MeV) 


Non relativistic 


12.254 


8.899 


7.724 


346 


Rclativistic 


11.334 


7.274 


4.56 


342 



TABLE V: The couplings G a , G u , G p of the non-relativistic (first line) and the relativistic (second line) version of the model 
for R B = 0.8fm and the m a = lOOMeV 



In a Fermi gas we have 

(2.TIM 1 



< K m (a) >= I dk^k 2 +M 2 (a), < K{a) >= ]T < tf m (a) > , (51) 



with kp(m) the Fermi level of the isospin species m. Using the solution for da we get, using standard techniques 



< H > , . 1 / dK 



1 ^ f kp{m) - - 1 9 — — — 9 



2r7i 2 \ 9fi 

1 \ rk F (m) ^ ^ 

+ 77i mz / / dk\dk 2 - _ 

(2ir)S^J ( kl -k 2 ) 2 +Al da 



rK F ym) ^ ^ I Q y fl y 

/ ^^TT r , 2 ~ 2 W fc ? + M\a) — Jk\ + M\a). (52) 

Jo (fei - fc?) 2 + m 2 ocr v do- v 



Finally we must add the contributions < >,< V p > associated with u> and p exchange. However, since these 
interactions are purely 2-body there is nothing new with respect to our previous calculation. One finds 

= -f (ft, - P«) 2 - G p J] C W — ^ / dfcr / dk 2 o -• (54) 

The energy density (|52I53I54[1 and its non-relativistic version, given in Section Hvl correspond to the same model. 
However, as in any theory, the parameters depend on the approximation scheme and therefore we need to readjust the 
couplings GvjGujGp to reproduce the saturation properties for normal nuclear matter. We do this for Rb — 0.8fm 
and m a — 7Q0MeV. As can be seen in Table [Vj the couplings we get are significantly different from those obtained 
with the non-relativistic expansion. The most important effect arises from the non-relativistic approximation (|21|l to 
the self consistentency equation Q48[l for a. If we were to use the same coupling constants in both equations, we would 
find a value of a which, at the saturation point, differs by 15%. 

In Fig. |3 we show the binding energy Eb, pressure P, incompressibility Kjq and sound velocity V s of symmetric 
nuclear matter predicted by our model (QMC) and the Skyrme (SkM*) interaction. We recall the definitions: 



*-i-«>-'w- x "<-H*&>- ,55) 

As we have restricted the flavor composition to protons and neutrons, the curves must not be taken too seriously 
for densities above 2 ~ 3 times po, because there the appearance of hyperons is expected to become important. 
However, the comparison between QMC and SkM* is still meaningful. One sees that the behavior of the two models 
is qualitatively similar in the range of densities displayed in the plots. Of course, at much higher density the SkM* 
model would violate the causality limit (V s — 1), because of its non-relativistic nature. By contrast the QMC model, 
because of its relativistic formulation, gives V s < 1 at any density. The QMC binding energy curve is stiffer than 
that for SkM*. This is not unexpected in view of the high value of K^(po) given by QMC. However, the difference 
is not spectacular when one considers the full density range. This is even more obvious when one looks at the 
incompressibility curve as a function of density. In fact the curves cross somewhere above p — lfm~ 3 . 

The situation changes radically when one considers the case of neutron matter - as shown in Fig. 0] The QMC 
equation of state is considerably stiffer than that for SkM* above p ~ 0.2fm~ 3 . This difference is not too surprising 
because in the Skyrme interaction the isospin dependence appears only in the 2-body interaction through the xq 
parameter, while in the QMC model it is also generated by the non-linearity of the a field equation, as can be seen 
in Eg . 1)3611 |29j . Thus one may expect that the neutron stars predicted by the QMC model will have higher maximum 
masses than in the SkM* model. However one cannot draw firm conclusions until the effect of hyperons has been 
calculated. 
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E B (MeV) P(MeV fm" 3 ) 




FIG. 3: Comparison of the QMC and SkM* model for symmetric matter. 
VIII. CONCLUSION 

Starting from the QMC model, we have derived a density dependent effective nucleon-nucleon force which is not 
limited to low density. The application of this effective interaction to the properties of doubly magic nuclei and then 
to nuclei far from stability yields results in quite acceptable agreement with popular effective forces which are far 
more phenomenological and which have rather more parameters. This provides one with confidence in the underlying 
physical concepts of the QMC model, most notably the scalar polarizability of a hadron which results from the self- 
consistent response of the confined valence quarks to the mean scalar field in the medium. An important corollary 
is that one can immediately apply the model to derive appropriate effective hyperon-nucleon effective forces as a 
function of density with no additional parameters. Given the lack of experimental guidance on anything other than 
the A — N force at relatively low density and the apparent importance of hyperons in dense matter, this will be our 
next application of the methods developed here. We remain mindful that our present work has been built upon the 
MIT bag model, which involves a relatively crude characterization of confinement. To the extent that the effect of 
internal structure is controlled by the scalar polarizability this is unimportant. Nevertheless, it will also be interesting 
to see the effect of using more sophisticated quark models. 

To put our results in perspective, we have applied the present approach to a relativistic Hartree-Fock treatment of 
dense matter consisting of nucleons. In the case of symmetric matter we find an equation of state which is not very 
different from the one predicted by the Skm* force but in the case of neutron matter it is definitively stiffer. However 
to infer the implications for the neutron stars properties we need to investigate the changes induced when hyperons are 
allowed. Finally, given that we are starting with a quark model for the nuclcon itself, it will be extremely interesting 
find the point at which a phase transition to quark matter becomes favourable. In this work we have formulated the 
QMC model so that it can be used formally at any density but it is clear that it must break at some point. Whether 
the central density in the massive neutron stars is high enough to reach this point is a tantalizing question. 
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FIG. 4: Comparison of the QMC and SkM* model for neutron matter. 
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[28] Spin and flavor labels are understood. 

[29] Note that this equation can only be used at moderate density since it is a non-relativistic approximation. 



